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1— I Abstract: We calculate contributions to the photon and gluon magnetic dipole oper- 

^ ators that mediate 6 — )• ,37 and b ^ d'j transitions in the Randall- Sundrum model of a 

i-Q warped extra dimension with anarchic bulk fermions and a brane localized Higgs. Unlike 

I the Standard Model, there are large contributions to the left-handed b quark decays, pa- 

fl^^ rameterized by the Wilson coefficient Cj, due to the pattern of bulk fermion localization, 

i-S^ and sizable contributions from the gluonic penguins, Cg , through renormalization group 
mixing. Further, unlike the Randall- Sundrum result for ^ — )• 67, the unprimed Wilson 

^^ coefficients receive non-negligible contributions from the misalignment of the bulk fermion 

<^ spectrum with the Standard Model flavor sector. We compare the size of effects and the 

/^J constraints imposed by the branching ratios Bt{B — )■ Xsj) and (Br(i? — t- Xdj)) within the 

\^ minimal and the custodial model. Within the custodial framework, we study the effect 

Cf^ on a number of benchmark observables and find that Br(i3 — )■ Xg^i^ ^~) and the forward- 

^~; backward asymmetry in B ^ K* ji^ ii~ remain close to their Standard Model predictions. 

T— I On the other hand, there can be large enhancements of the time-dependent CP asymmetry 

^ in i? — 7- 7^*7 and the transverse asymmetry Aj.' . 
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1 Introduction 

The Randall-Sundrum (RS) scenario of a warped extra dimension provides an elegant so- 
lution to the hierarchy problem [1-6] and a way to understand strongly coupled dynamics 
through the AdS/CFT correspondence [7-9]. For reviews see [10-13]. One of the promising 
phenomenological features to come out of this framework is an explanation of the Stan- 
dard Model (SM) flavor structure through the split-fermion scenario [3, 5, 14, 15]. In these 
models the Yukawa matrices are anarchic and the spectrum of fermion masses is gener- 
ated by the exponential suppression of zero mode wavefunctions with a brane-localized 
Higgs [16]. This also automatically generates hierarchical mixing angles [16-18] and sup- 
presses many tree-level flavor-changing neutral currents (FCNCs) through the RS-GIM 
mechanism [16] . In order to protect against large contributions to the T parameter coming 
from bulk gauge fields, one may introduce a gauged custodial symmetry [19] that is broken 
on the boundaries; a straightforward discrete extension of such a symmetry also protects 
against corrections to the Zbh vertex [20, 21] and flavor changing couplings of the Z boson 
to left-handed down-type quarks [22, 23]. 

These flavor protection mechanisms are not always sufficient to completely protect RS 
models from stringent experimental flavor constraints. In the quark sector, the tree-level 
exchange of Kaluza-Klein (KK) gluons and neutral electroweak gauge bosons contributes to 
meson- antimeson mixing and induces left-right operators. These operators are not present 
in the SM and receive a signiflcant enhancement through QCD effects due to their large 
anomalous dimension. In the kaon system they are also chirally enhanced by a factor of 
m^/rrig. These contributions lead to new CP violating effects in the kaon system, namely 
the well-measured observable ek, and result in generic bounds of 0(10 — 20 TeV) for the 
KK gluon mass [23-29]. To reduce this bound, one must introduce additional structure 
such as horizontal symmetries [30, 31], flavor alignment [32, 33], or an extended strong 
sector [34]. Alternately, one may promote the Higgs to a bulk fleld [35] to localize the 
fermion zero modes closer to the UV brane. 

Additional constraints on the RS flavor sector come from loop-induced dipole operators 
through penguin diagrams. The flrst estimates for these operators were performed in [16, 
17, 36, 37] assuming UV sensitivity at all loops within the 5D effective theory and a 
calculation within the two-site approach was performed in [35]. In [16] the bound Mkk > 



0(10 TeV) was derived from the constraint on the neutron electric dipole moment. The RS 
dipole contributions lead to dangerously large effects in direct CP violation in the K — )• vrvr 
decays measured by the ratio e'/e [38]. Combining the bound from the latter ratio with the 
£k constraint leads to a lower bound on the KK scale independent of the strength of the 5D 
Yukawa. More recently it was shown that even for the brane Higgs scenario the one-loop 
induced magnetic penguin diagrams are finite in RS and can be calculated effectively in a 
manifestly 5D formalism [39]. The lepton flavor violating penguin /i — )• 67 sets bounds on 
the KK and anarchic Yukawa scales that are complementary to tree-level processes, so the 
tension between these bounds quantifies the degree of tuning required in the 5D Yukawa 
matrix [37]. 

In this paper we examine the calculation and phenomenological observables of the 
quark sector processes 6 — )• (77 (g = s, d) in the RS framework with a brane-localized Higgs 
field using the mixed position-momentum space formalism. These processes differ from 
their leptonic analogs for various reasons beyond the spectrum and diagrams involved. 
Firstly, while the branching ratio of /U — )• 67 is only bounded from above, the branching 
ratios for B — )• Xg^ and, to a lesser extent, B — )• X^j are well- measured and in good 
agreement with the SM. Secondly, theoretical predictions are more involved due to the 
renormalization group (RG) evolution from the KK scale to the B meson scale and hadronic 
effects at the latter scale. The RG running over this large range of energy scales introduces 
a sizable mixing between the various effective operators, so that one must also include the 
effects of the magnetic gluon penguin Cg in addition to the magnetic photon penguin C7 . 

After reviewing the flavor structure of RS models in Section 2, we calculate the Cj 
and Cg Wilson coefficients of the quark dipole operators in Section 3. We provide explicit 
formulae for the dominant RS contributions to the Wilson coefficients at the KK scale in 
both the minimal and custodial models and analyze the size of these contributions. In 
Sections 4 and 5, we subsequently perform the RG evolution down to the B meson scale 
and obtain predictions for the branching ratios Br(i? — )• Xg^dj)- 

Finally, in Section 6, we investigate the phenomenological implications on a number 
of benchmark observables related to the photon and gluon penguin operators. We first 
show that these operators give non-negligible constraints for both minimal and custodial 
models. We then restrict our attention to realistic models with a bulk custodial symmetry 
SU{2)l X SU{2)ji X U{l)x X P^/j and consider the effect of benchmark observables on 
points in parameter space that pass tree-level constraints as evaluated in [25] . Rather than 
performing a detailed analysis of all observables provided by the B — )• X^'j, B — )• K*j, 
B — )• XsH'^ iJL~ and B — )• K*fj,~^fj,~ decay modes, we focus on a number of benchmark 
observables in order to illustrate the pattern of effects and leave a more detailed analysis 
for future work. Specifically we study: 



• 



The branching ratio Br(i? — t- Xg^) and the CP averaged branching ratio (Br{B — )• 
Xfi't)) which we impose as constraints on our parameter scan. 

The branching ratio Br(i? — t- Xsfi^fi^) and the forward backward asymmetry ^fb 
in i? — )• K* 11^ iji~ . Stringent data that are in good agreement with the SM exist for 
both observables, placing strong bounds on various new physics (NP) scenarios. The 



custodial RS model naturally predicts small effects in these observables since they are 
rather insensitive to NP contributions to the primed magnetic Wilson coefficients. 

• The time-dependent CP asymmetry Sk*--! in B — )• K*^ and the transverse asymmetry 
Arp in i? — )• K*fj,^fi~, evaluated in the region of low dimuon invariant mass 1 GeV'^ < 
q^ <6GeV2. 

Since the RS contributions generally exhibit the hierarchy AC'^ ^ AC7 [16, 36] the 
latter observables are particularly suited to look for RS contributions. CP asymmetries in 
radiative B decays were already suggested in [16, 36] as good probes to look for RS effects. 
We quantify the possible size of effects and study the possible RS contributions to the 
various observables in a correlated manner. We also included the transverse asymmetry 
A)p \ which has not been considered in the context of RS models before. 

2 Flavor in Randall- Sundrum models 

We summarize here the relevant aspects of flavor physics and the RS scenario. For a 
review of the general framework see e.g. [11-13, 23, 28]. We consider a 5D warped interval 
z € [R, R'] with an infrared (IR) brane at z = i?' ~ (TeV)~^ and an ultraviolet (UV) brane 
at z = -R ~ Mpi, the AdS curvature scale. In conformal coordinates the metric is 

ds'^ = (-\ {dXf.dx^T]''^ - dz^). (2.1) 

One may recover the classic RS conventions with the identifications z = Rexp{ky) and 
k = 1/R, kexp{—kL) = 1/R'. Fermions are Dirac fields that propagate in the bulk and 
can be written in terms of left- and right-handed Weyl spinors x &nd ■i/' via 

In order to obtain a spectrum with chiral zero modes, fermions must have chiral (orbifold) 
boundary conditions, 

M^^',R)=M^^',R') = and xr{x^, R) = Xr{oo'',R') = 0, (2.3) 

where the subscripts L and R denote the SU(2)l doublet (L) and singlet (R) representa- 
tions, i.e. the chirality of the zero mode (SM fermion). The localization of the normalized 
zero mode profile is controlled by the dimensionless parameter c, 

x(°)(x,z) = -^(|)'(|^)"Vcxi°Hx) and ^^(x, z) = xS(x, .), (2.4) 



^\RJ \R'^ 

where c/R is the fermion bulk mass. Here we have defined the RS flavor function charac- 
terizing the fermion profile on the IR brane. 



We assume that the Higgs is locahzed on the IR brane. The Yukawa couphng is 

1 



Syuk = I (fx [ — 



^ 



Qi-HR Yu,ijUj + Qi-HR Ya^ijDj + h.c, 



(2.6) 



where Yij are dimensionless 3x3 matrices such that (Xb)ij = RYij is the dimensionful 
parameter appearing in the 5D Lagrangian with Y assumed to be a random 'anarchic' 
matrix with average elements of order Y^,. After including warp factors and canonically 
normalizing fields, the effective 4D Yukawa and zero mode mass matrices are 



yfl" = fc,Y,,f.,,^ m.^ = —yff^, (2.7) 



so that the fermion mass hierarchy is set by the /i ^ /2 ^ fs structure for both left- and 
right-handed zero modes. At the same time, the hierarchical pattern of the CKM matrix is 
also generated naturally. In other words, the choice of c for each fermion family introduces 
additional flavor structure into the theory that generates the zero mode spectrum while 
allowing the fundamental Yukawa parameters to be anarchic. 

In this document we work in the gauge basis where the bulk mass matrices and the 
interactions of the neutral gauge bosons are flavor diagonal but not flavor universal. The 
Yukawa couplings are non-diagonal in this basis and cause the resulting fermion mass matri- 
ces to be non-diagonal. Since these off-diagonal entries are governed by the small parameter 
vR' , we will treat them as a perturbative correction in the mass insertion approximation. 

Realistic RS models typically require a mechanism to suppress generically large con- 
tributions to the Peskin-Takeuchi T parameter and the Zbb coupling; a common technique 
is to extend the bulk gauge symmetry to [19-21, 40-43] 

SU{3)c X SU{2)l X SU{2)n x U{l)x x Plr. (2.8) 

Here Plr is a discrete symmetry exchanging the SU{2)l and SU{2)ji factors; in order to 
protect the left-handed Zbb coupling from anomalously large corrections, the left-handed 
down type quarks have to be eigenstates under Plr- This in turn requires enlarged fermion 
representations with respect to the minimal model. Specifically the quark representations 
containing the SM zero modes are {i = 1,2,3): 



xl'{-+h/3 qT (.++)2/3 
,xJ(-+)2/3 qii++)^i/3, 



ci _ Al I— r;5/3 Vl V-r-r;2/3 ,^ qx 



-,2R 



->3R 



2/3 
nM++)2/3, (2.10) 

/ ^U-+h/3 \ ( <(- + )5/3 \ 

nR®Tij,= \ u'j^{-+)2/3 e U'^K-+)2/3 ■ (2.11) 

Here ^]^^ is an SU{2)l x SU{2)r bidoublet, S,2r is singlet under both SU{2)s, and Ti^j^ and 
TIj^ are triplets under SU{2)l and SU{2)r respectively, with all of them carrying U{l)x 
charge +2/3. The corresponding states of opposite chirality are obtained by reversing the 
boundary conditions. 



As we will see later, while the additional gauge bosons present in the custodial model 
do not have a significant impact on the 5 — >• 57 and b ^ qg {q = d, s) amplitudes, the 
additional fermion modes contribute and generally enhance the effect. 

3 Calculation of the b — )■ qj Penguin in RS 

We now calculate the RS contributions to the b ^ qj and b ^ qg (q = d, s) decays. These 
contributions are calculated at the KK scale Mkk ~ 1/-R'; in subsequent sections we will 
relate these to renormalization group (RG) evolved coefficients and observables at the low 
scale ~ nib. 

We only evaluate the dominant diagrams, working in Feynman gauge and the mass 
insertion approximation, where the expansion parameter is vR'/^/2 ~ 0(0.1). We have 
checked explicitly that the diagrams presented here dominate those that were neglected 
by at least an order of magnitude; a more detailed calculation is beyond the scope of this 
work and, in our opinion, premature before the discovery of RS KK modes. We refer 
to [39] for details of the 5D calculation, Feynman rules, and guidelines for estimating 
the dominant diagrams. For additional notation and conventions, especially with respect 
to the custodially protected model, see [23]. See Appendix D for comments on theory 
uncertainties. 

3.1 Effective Hamiltonian for 6 — )■ 97 transitions 

The 6 — )• (77 (g = d, s) transitions are most conveniently described by an effective Hamilto- 
nian in the operator product expansion, see e. g. [44] for a review. The dipole terms most 
sensitive to new physics are 



V2 



'HeS= — ^Vt^Vtb 



C7{f^)Q7if^) + C'^{^i)Q',{^,) + Cs{fs)Qs{fi) + C'sifi)Q's{f^)\ + h.c, (3.1) 
where we neglect terms proportional to Vuq^ub- The effective operators are 

Q7 = £^mb{qa^,PRb)F^''' Q', = ^rrn,{qa^,PLb)F>''' (3.2) 

^8 = {^rnb{qa^,T^PRb)G^'^^'' Q's = |^mfe(ga^,T'^PL6)G'''^'^ (3.3) 

where Pl,r = (l-F7^)/2. In this document we will focus on new contributions from the RS 
model to these operators. There are also contributions from non-dipole operators Qi^...fi 
and their chirality-flipped (primed) counterparts, but these are far less sensitive to NP and 
can be assumed to be equal to their SM contributions^. 

At leading order in the SM, the primed Wilson coefficients C^ g are suppressed by 
ms/nib and therefore negligible, so the relevant Wilson coefficients at the scale M^y are 

C^'^iMw) = -\d',{x,) , Cf^{Mw) = -lE',{xt), (3.4) 



^The impact of flavor changing neutral gauge bosons on the operators Qi....,6 has recently been studied 
in [45]. Since the relevant contributions in RS are suppressed both by the KK scale and the RS GIM 
mechanism, the contributions are expected to be small and will be neglected in this paper. 



where xt = rn^/Myy, and Do(^i) ^ 0-37 and -E'o(^<) ^ 0.19 are loop functions given 
explicitly in (3.15-3.16) of [46]. In what follows we refer to the RS contributions to these 
operators as ACj § . 

3.2 Structure of the amplitude 

In order to calculate the b — )• (s, d)j and b — )• (s, d)g penguins, we work in a manifestly 5D 
framework. Unlike the 4D KK reduction, this procedure automatically incorporates the 
entire KK tower^ at the cost of an expansion with respect to the Higgs-induced mass term 
{r^vR'). 

Using the on-shell condition for the photon, the general form of the left-to-right chi- 
rality //"(p) — )• ffip')"/ amplitude, C7, in a 5D theory can be written as [39, 48] 

^L^-^B, = ^ ^E ( «^-^4^fc^4 + ^^^■4) fQjD.uP [{p + p'r - {m, + m,)r] u^,e, 

(3.5) 

where e is the photon polarization. The chirality flipped amplitude is given by the conjugate 
of this result, Aij^i^ij = (A^/^j_^/j»)^. The expression for the gluon penguin is analogous 
with the appropriate substitutions. Using the fermion equations of motion, the term in the 
square brackets gives the required dipole structure a^^F^^-, so a simple way to identify the 
gauge-invariant contribution to the amplitude is to determine the coefficient of the {p+p')^ 
term [48] . In [39] this observation was used to show the manifest one- loop finiteness of these 
dipole transitions in 5D theories. Matching (3.5) to the effective Hamiltonian (3.1) yields 
expressions for the RS contributions to the Wilson coefficients, AC. 

We refer to the coefficients a^i and bij in (3.5) as the anarchic and the misalignment 
contributions, respectively. They are products of couplings and dimensionless integrals 
whose flavor indices reflect the bulk mass dependence of internal propagators. Upon di- 
agonalizing the SM fermion mass matrix, the anarchic term a is not diagonalized and 
generally remains anarchic. On the other hand, in the limit where the bulk masses are 
degenerate, the flavor structure of the b term is aligned with the SM Yukawa matrices and 
thus contains no flavor-changing transitions in the mass basis [16, 35, 37]. This alignment 
is pronounced for the first and second generation fermions because their bulk masses are 
nearly degenerate, but special care is required for the third generation quarks since these 
are localized towards the IR brane. The physical contribution of the b coefficient comes 
from the robustness of off-diagonal elements of bijYijfqJo. after passing to the basis in 
which YijfQ^fj:). is diagonalized. Contrary to the usual assumption of Yukawa anarchy, 
the overall size of the b term depends on the misalignment of the specific anarchic Yukawa 
matrix relative to the set of bulk masses as flavor spurious. One expects diagrams with 
internal zero modes to give the dominant contributions, since these are the most sensitive 
to the bulk mass spectrum and hence robust against diagonalization; this intuition is con- 
firmed by our numerical scans. One measure of this effect is the la standard deviation from 



^ An alternate method of including the entire KK tower based on residue theorems was presented in 
[47], though it obfuscates the physical intuition presented below. 



6 = in a scan over random anarchic matrices [39]; we use this to identify the dominant 
contributions to this misahgnment term. 

By assumption, the anarchic contribution is independent of the SM flavor sector, so 
there is no analogous alignment suppression to the a coefficient. However, depending on 
the internal modes in the loop, each diagram contributing to this term carries one of two 
possible independent flavor spurious that can be built out of the Yukawa matrices that 
may enter this product: YuYuYJ and YjVdYj. These matrices may have arbitrary relative 
phase, so the two terms may add either constructively or destructively. The misalign- 
ment contribution is a third independent flavor spurion, which also carries a relative phase 
dependent on the particular choice of parameters. 

We express the anarchic (a) and misalignment (6) coefficients in terms of dimension- 
less integrals, which are defined in Appendix A. To explicitly demonstrate the calculation 
of diagrams in the 5D mixed position/momentum space formalism, we present a sample 
calculation of the anarchic contribution to Cj in Appendix B. The Cg diagrams where a 
gluon is emitted from an internal gluon have integral results that are typically 0{1) while 
the integrals for the other diagrams are typically 0(10"^) in magnitude. Note that the 
contribution to a from each diagram matches what is expected from a naive dimensional 
analysis. This is in contrast to the analogous calculation for /i — )• 67, where the leading 
diagrams are smaller than the naive estimated size. There are thus no problems with the 
two-loop contribution yielding a larger contribution than expected from the perturbative 
expansion. 

Below we present the calculation for the right-to-left chirality (unprimed) Wilson coef- 
ficients ACy^g for 6 —7- (7; the left-to-right chirality (primed) Wilson coefficients are obtained 
by Hermitian conjugation of the g — t- 6 amplitude. The anarchic contribution to the left- 
to-right chirality coefficients are enhanced over the right-to-left coefficients by a factor of 
fbi^/fbii: while the misalignment contribution is of the same order of magnitude. This be- 
havior is explained qualitatively in Appendix C and demonstrated numerically in Section 6. 

3.3 Calculation of Acj'^ 

Fig. 1 shows the dominant contributions to the Cj photon penguin operator. The RS 
contribution to the b ^ q'j Wilson coefficient is 



AC7 



-vR 



/2 



Sm^Gp 



(VtlVtb) 



ijki 



qi I 
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D, 



J2a,eYr,%-,Y,f + b,,Y, 



ij -' ij 



kl 



U, 



jb 



(3.6) 



U '^•^ are the rotation matrices between the 5D gauge and the light down quark mass 
bases. 

Note that throughout our analysis we use the tree level matching condition for the 
5D gauge couplings and neglect possible brane kinetic terms that may alter this matching. 



(') 



.(') 



While this affects the misalignment contribution to C) and the calculation for Cg , the 
anarchic contribution to C\ , containing only one gauge coupling vertex instead of three, 
remains relatively unaffected. Since the latter gives the dominant contribution to the 





H- 

(a) Charged Goldstone loop 



G G 

(b) Gluon (G,! or G5) loops with a single mass insertion 



Figure 1. Leading contributions to the anarchic (a) and misalignment (6) terms of the C7 Wilson 
coefficient. Arrows indicate SU(2)l representation; this is equivalent to labeling the chirality of 
the zero mode for SM fields. Here Q, U and D denote the 5D chiral fermion fields containing 
the SM left-handed doublets and right-handed up and down singlets, respectively. H^ is the 
charged component of the Higgs doublet that serves as the Goldstone boson of W^ after electroweak 
symmetry breaking, and G is the 5D gluon field. Additional diagrams related by exchanging the 
order of the mass insertion and photon emission are left implicit. 



observables discussed in section 6, we do not expect this assumption to have a significant 
impact on our predictions. 

3.3.1 AC7: anarchic contribution 

The dominant anarchic contribution is the diagram with one mass insertion and a charged 
Higgs (Goldstone) in the loop, Fig. la. Note that this diagram is not present in the 
analogous leptonic penguin, which has a neutrino in the loop. The analogous diagram 
with the photon emitted off the charged Higgs propagator is found to be suppressed by a 
factor of {m]yR')'^ ~ 10~^ due to an algebraic cancellation [39], while the one with the mass 
insertion on an external fermion leg is suppressed by ruq/v since the external brane-to-brane 
fermion propagator must be a zero-mode. All other diagrams contain two additional mass 
insertions — necessary to obtain the required structure of a product of three Yukawas — 
and are therefore also suppressed by a factor of {vR' /V2)^ ~ 10~^. Of these neglected 
diagrams, the next-to- leading diagrams contributing to this coefficient are gluon loops with 
three mass insertions, which carry a gauge coupling enhancement of glhiB! /R w 36 but 
are suppressed due to the two additional mass insertions, the quark charge {Qd = —1/3), 
and different topologies; they are ~ 5% corrections to the leading contribution. Note that 
these diagrams carry an independent flavor structure (l^'Y^Yj) and can interfere either 
constructively or destructively with Fig. la. 

The value for the a coefficient in (3.6) coming from the penguin in Fig. la is a dimen- 
sionless integral whose explicit form is given in (A. 11), 



a = QuICva J 

where Qu = 2/3 is the charge of the internal up-type quark. 



(3.7) 



3.3.2 ACt: misalignment contribution 

The dominant misalignment contributions come from gluon diagrams with a single mass 
insertion. As shown in Fig. lb, this insertion can either be on an internal or external fermion 
line. All other diagrams contain electroweak couplings and hence are subdominant. The 
final misalignment contribution in (3.6) is 



Qa 



''^""ll 



Ic,, 



(3.8) 



Here Qd is the charge of the internal down- type quark, 4/3 is a color factor, In R'/R is 
a warp factor associated with bulk gauge couplings, and Iq^^ is a dimensionless integral 
defined in (A. 12). 

3.4 Calculation of AC^'^ 

The gluon penguin operators Cs and Cg differ from their photon counterparts due to 
additional QCD vertices available and the magnitude of the QCD coupling, g^j-)/R = 
g^ In R'/R ~ 36. Because of this, the dominant diagrams contributing to b ^- qg cannot 
be obtained from 6 — )• 57 by simply replacing the photon with a gluon in the leading 
diagrams for Cj . The general expression for AG's is the same as that for AC7 in (3.6), 
with coefficients a and b coming from the diagrams shown in Fig. 2. 




(a) Charged Goldstone and three mass insertion gluon loops. Not 
shown: gluon loop with two and three internal line mass insertions. 




(b) One mass insertion gluon 
loop 



Figure 2. Leading contributions to the a and b terms of the Cs Wilson coefficient following the 
notation of Fig. 1. G^ refers to only the gluon four- vector. 



3.4.1 ACs: anarchic contribution 



.(') 



There are two classes of dominant contributions to the anarchic (a) coefficient in Cg . In 
addition to the charged Higgs diagrams analogous to Fig. la, there are gluon diagrams 
with three mass insertions on the fermion lines, which are now sizable due to the size of 
the strong coupling constant and the three-point gauge boson vertex (as mentioned earlier, 
the dimensionless integral associated with this digram is 0(1), while all other diagrams 
have 0(0.1) integrals). Of the latter class, one only needs to consider diagrams with at 
most one mass insertion on each external leg since sequential insertions on an external 



leg are suppressed by factors of rUqE! . Note that these two sets of diagrams contribute 
with different products of Yukawa matrices; while the Higgs diagrams are proportional to 
YuYuYJ, the gluon diagrams are proportional to YjYdYj. Thus these two terms may add 
either constructively or destructively and may even add with different relative sizes if there 
is a hierarchy between the overall scale of the up- and down-type 5D anarchic Yukawas. 
The a coefficient is 

where we have written © to indicate that the two terms carry independent flavor spurious. 
Here Ic^a is the same dimensionless integral appearing in (3.7). The second term includes 
color factors, warped bulk gauge couplings, and explicit mass insertions in addition to the 
dimensionless integral Icg^ defined in (A. 15). 

3.4.2 ACg: misalignment contribution 

The single mass insertion gluon emission diagram in Fig. 2b gives the dominant misalign- 
ment term. Additional diagrams with the gluon emission from the quark line are suppressed 
by a relative color factor of 1/6 versus 3/2 and can be neglected. Diagrams with a scalar 
(G5) gluon or the mass insertion on an external leg do not carry an internal fermion zero 
mode and become negligible after rotation to the mass basis, as discussed earlier. Diagrams 
with electroweak gauge bosons in the loop are suppressed due to the smaller size of the 
gauge coupling. The expression for the dominant diagram is 



5'ln^ Ucs,- (3.10) 



with Icg^ defined in (A. 19). We have again pulled out an explicit color factor and the 
warped bulk gauge coupling. 

3.5 Modifications from custodial symmetry 

In models with a gauged bulk custodial symmetry, the additional matter content may 
also contribute to the b — )• q^{g) transitions. By construction, boundary conditions for 
custodial fermions are chosen such that they have no zero modes. The misalignment (b) 
coefficients do not receive any significant corrections from custodial diagrams: diagrams 
with custodial gauge bosons are suppressed due to electroweak couplings, while those with 
custodial fermions do not carry internal fermion zero modes and become negligible after 
rotation to the mass basis. 

The leading custodial contributions to the anarchic (a) coefficients are shown in Fig. 3; 
these are the same diagrams that contribute to the anarchic (a) terms of the C7 and Cg 
Wilson coefficients and now appear with additional custodial fermions, denoted by U\ 
U" , and D' . These are the only custodial diagrams that give contributions comparable to 
those in Fig. 1 and Fig. 2. The remaining diagrams consist of W and Z loops, which, as 
mentioned earlier, are suppressed by a factor of ~ 10^^ relative to the Higgs loops due 
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Figure 3. Additional custodial diagrams contributing to the C7 and Cs coefficients. 



to the two additional mass insertions, and remain negligible despite the larger multiplicity 
due to the extended electroweak sector. 

Since the custodial fermions C/', C/", and D' have the same IR boundary condition as 
their SM counterparts but the opposite UV boundary condition, and since the localization 
of the Higgs pulls the loop towards the IR brane, the contribution of these custodial 
diagrams is well- approximated by the contributions of their SM counterparts. Since the 
minimal model diagrams are dominated by the KK fermion contribution, it is reasonable 
that the custodial modes should contribute approximately equally to the process. 

Observe that each of these custodial contributions is proportional to YjYdVj. In 
particular, the custodial Higgs diagrams carry a flavor structure that is independent of 
that of their minimal model counterparts. Also, note that the U' and U" couplings to 
the charged Higgs come with a factor of l/-v/2 while the D' coupling to the Higgs does 
not [25]. Thus the additional custodial diagrams contribute an analytic structure that is 
nearly identical to the minimal model diagrams except for the Yukawa matrices, which now 
come with the product Y^YdYj. Since this is independent of the Y^YuYu flavor spurion in 
the minimal model diagrams, the addition of the custodial diagrams generically enhances 
the penguin amplitude by less than the factor of two that one would obtain in the limit 
Yd = Yu. This shows that while custodial symmetry can be used to suppress tree- level flavor 
changing effects in RS models, this comes at the cost of generically enhancing loop- level 
flavor processes. 



4 Radiative B decays 

We now examine the physical observables most directly related to the parton-level b 
5(7, g) operators derived above: B meson decays with an on-shell photon. 

4.1 The B -)■ Xs,dl decay 

The SM predictions for the inclusive decays B — )• X^^dl are [49, 50] 



+2.6\ 

-3.1^ 



10" 



Br(B ^ X,7)sM = (3.15 ± 0.23) • 10"^ , (Br(B ^ Xdl))sM = (15.4 

These can be compared to the measured values [51] 
Br(S ^ X,7)exp = (3.55 ± 0.27) • 10"^ , (Br(5 ^ X^^)),^^ = (14 ± 5) • lO''^ 



(4.1) 
(4.2) 
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Here {Br{B — t- Xdj)) refers to the CP averaged branching ratio in which the hadronic 
uncertainties cancel to a large extent [52]. We have extrapolated the experimental value 
for (Br(i? — )■ X^'y)) to the photon energy cut E^ > 1.6 GeV used for the theory prediction. 
Rather than performing an extensive error analysis, we simply require the new RS 
contributions to fulfill the constraints 



ABt{B ^ Xs7) = 


Br(5^X,7)exp - 


- Br{B^Xs7)sM -- 


= (0.4±0.7)-10-\ (4.3) 


ABr(S -^ Xdj) = 


(Br(B ^ Xd7))cxp - 


- {Bi{B -^ Xd-f))sM ■ 


= -(l±ll)-10-^ (4.4) 



Neglecting all uncertainties associated with NP contributions, these constraints represent 
the 2(7 ranges when combining experimental and theoretical uncertainties in quadrature. 
Although the data and prediction for B — )■ X^j are currently less precise than those for 
B — )• ^^7, an important and partly complementary constraint can be obtained from the 
former decay, as recently pointed out in [50] . Since the data for B — )• X^j lie slightly below 
the SM prediction, ABr(B — )• X^j) < is somewhat favored, leaving little room for NP 
contributing to Cj. In contrast, a positive NP contribution to Br(B — )■ ^^7) is welcome to 
bring the theory prediction closer to the data. We note that if the tree level values for the 
CKM parameters are used instead of the SM best fit values, the predicted central value for 
{Bt{B — )• Xd7))sM rises to about 19 • 10^^, increasing the tension with the data. 

4.2 Master formula for Br{B -^ Xgf) 

Following the strategy of [45, 53, 54], which use the results of [55], the "master formula" 
for the inclusive B — )■ Xs^ branching ratio in terms of the SM branching ratio, BrgM, and 
NP contributions to the Wilson coefficients is 

Br(B ^ Xs-f) = BrsM + 0.00247[|AC7(/i6)|2 + \AC'j{fib)\'' - 0.706 Re(AC7(/i6))] . (4.5) 

The RS contributions to AC^ (fib) are obtained from the RG evolution of AC^ and ACg , 
calculated in Section 3 at the high scale Mkk = 2.5 TeV, down to the B scale, fib = 2.5 GeV, 

ACi'\fib) = 0.429 AC}'^ (Mkk) + 0.128 AC^'^ (Mkk) • (4.6) 

All known SM non-perturbative contributions have been taken into account while the RS 
contribution is included at leading order neglecting uncertainties. This approach is an 
approximation to studying the effects of RS physics on the decay in question; however, in 
view of the other uncertainties involved — such as the the mass insertion approximation and 
taking into account only the leading diagrams — this approach gives sufficiently accurate 
results to estimate the size of RS contributions. A more accurate and detailed analysis is 
beyond the scope of our analysis and, in our view, premature before the discovery of RS 
KK modes. 
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4.3 Master formula for (Br(S -)• Xdi)) 

A master formula can be obtained in a similar manner for the CP-averaged B — ?• X^^ 
branching ratio. Using the expressions collected in [50, 54, 56] we find 

(Br(B ^ Xd-i)) = (BrsM> + 10"^ [l.69 (| AC7P + |AC^|2) + 0.24 (| ACgp + jAC^p) 

+ 1.06 Re [ACyAC^ + AC^AC^*] - 3.24Re(AC7) 
-0.16Im(AC77) - 1.03 Re(AC78) -0.04 Im(AC8)j (4.7) 

where all of the RS contributions to the 6 — t- d Wilson coefficients /S.C\ g are evaluated at 
Mkk- 

4.4 Analytic estimate of constraints 

Assuming anarchic Yukawa couplings, one may estimate the size of the RS contributions 
to the Wilson coefficients in terms of the anarchic coefficients in Section 3.3.1, 

|AC7(Mkk)'^^''^^| ~ 1^/^*^' ~ 0.015 an^ (h^) ' ^^-^^ 

|AC^(Mkk)^-^^| ~ -^aY^H^^^^ ~ 0.18 an^ f — ^)' , (4.9) 

|AC^(Mkk)'^''^| ^aY^E!^ ""f „ ~ 0.20 a^^ ( ^' \ , (4.10) 

where we neglect the misalignment contributions. Here 1^ is the average size of the anarchic 
Yukawa couplings Yij which we assume to be equal for Yu and Y^. 

Generically the contribution to the chirality-ffipped operator Cy is larger than the one 
to C7 by more than an order of magnitude. This is a direct consequence of the hierarchical 
pattern of quark masses and CKM angles: in order to fit the observed spectrum, the left- 
handed bi quark has to be localized close to the IR brane, and consequently its flavor 
violating interactions are far more pronounced than those of the right-handed bR. 

Neglecting the subdominant contributions from AC7 and ACg , we can constrain the 
size of AC7 by making use of the data on Bic{B — )• Xs^) and {Bic{B — )• Xdj))- We obtain 
the following constraints from the master formulas and the experimental constraints quoted 
above: 

|AC;(Mkk)*^'^| <0.47, |AC:^(Mkk)''^'^^| <0.77. (4.11) 

Using (4.9-4.10) and a ~ 0.33 we can derive an upper bound on the size of the Yukawa 
couplings, y*, 

— ^ — J < 2.8 from B -^ Xsi , (4.12) 

—^ — - < 3.4 from B -^ Xaj , (4.13) 

For R' = 1 TeV^ these are of the same order as the perturbativity bound on the Yukawa 
coupling [24]. We see that the generic constraint from B — ?■ Xgj is slightly stronger than 
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that from B — )■ Xd'j due to the larger uncertainties in the latter case. However, since they 
only differ by an C(l) factor, in specific cases the latter constraint may be more restrictive, 
so one must take both processes into account when constraining the RS parameter space. 

4.5 CP asymmetry in B —^ K*'~f 

Like many extensions of the SM, RS generally induces large CP violating phases. It is 
thus of great interest to also study CP violation in 6 —t- 57 transitions. While the direct 
CP asymmetry in the inclusive B — )■ X^^ decay is in principle highly sensitive to NP 
contributions, in practice the SM contribution is dominated by long-distance physics and 
therefore plagued by large non-perturbative uncertainties [57]. Consequently, a reliable 
prediction in the presence of NP is difficult. 

Fortunately, a theoretically much cleaner observable is provided by the B — t- K*^ 
decay. While its branching ratio is plagued by the theoretical uncertainty of the B — t- K* 
form factors, this form factor dependence largely drops out of the time-dependent CP 
asymmetry [58-60] 

;-, ^-^^ y-^. w^ = SK*-y sin AMrfi) - C;^*^ cosfAM^t . (4.14 

The coefficient Sx'-f is highly sensitive to new RS contributions. At leading order it is 
given by [59, 61] 

S^'^ ^ ^q^lr^rplm (e-^^'^CrC;) , (4.15) 

where the Wilson coefficients are to be taken at the scale ^t. (f)^ is the phase of B^~B^ 
mixing, which has been well measured in B^ — t- J/ipKs decays to be sinc/)^ = 0.67 it 0.02 
[51]. 

From (4.15) we see that SK*'y is very sensitive to new phsyics in the chirality flipped op- 
erator C'y and vanishes in the limit Cj — )• 0. Consequently the SM prediction is suppressed 
by the ratio ms/nn, and is therefore very small [60], 

5|¥^ = (-2.3±1.6)%. (4.16) 

Measuring a sizable CP asymmetry S'x*7 would thus not only be a clear sign of physics 
beyond the SM, but unambiguously indicate the presence of new right handed currents. 
The present experimental constraint [51, 62, 63], 

5^-P^ = -16%±22%, (4.17) 

is still subject to large uncertainties but already puts strong constraints on NP in 6 — )• s 
transitions [61]. A signiflcant improvement is expected soon from LHCb, and the next 
generation B factories will reduce the uncertainty even further. 

5 Semileptonic B decays 

Semileptonic B decays such as i? — t- Xgix^ ii~ and B — )■ K* fi^ ijl^ offer an interesting 
opportunity to not only look for deviations from the SM, but also to identify the pattern 
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of NP contributions and therewith distinguish various NP scenarios. These decays receive 
contributions from semileptonic four-fermion operators {sb){fifj,) in addition to the magnetic 
dipole operators discussed earher. While the dipole operators receive RS contributions first 
at the one- loop level as required by gauge invariance, the four fermion operators are already 
affected at tree level by the exchange of the Z boson and the heavy electroweak KK gauge 
bosons. 

In this section we discuss the effective Hamiltonian for b — )• sfi~^^~ transitions. Subse- 
quently we will review a number of benchmark observables that are relevant for the study 
of RS contributions. 

5.1 Eflfective Hamiltonian for b —^ sfi~^fi~ transitions 

The effective Hamiltonian for b — )• s^^ ^^ reads 

T-Lcs = T-Lcsib -^ s-f) ■j^Vtl'^tb Cw{^^)Qw{^J') + Cgv(p)Q9v(.f^) 

-FCioA(^)QlOA(/^) + C'ioA(/i)<5lOA(M) +h.c. , (5.1) 

where we neglect the terms proportional to V*gVub, and 

Qgv = 2{s-f^PLb)ifljy) Q'gy = 2{s^^PRb){fl^f'fi) (5.2) 

Qioa = 2(s7a.^l6)(/x7'^75/^) QIqa = 2(s7^P/?6)(/27''75^)- (5-3) 

In the SM only the unprimed Wilson coefficients are relevant. At the scale Mw they are 
given by 



Cl^iMw) = ^ 



4:Zo(Xt) 



sin 6w 



ySMn^^ \ _ « Yo{xt) 



C^^-^(Mw) = --^^^ (5.4) 

zvr sm Ow 



where xt = ml/My^ and the dimensionless loop functions YQ{xt) ~ 0.94 and ZQ{xt) ~ 0.65 
are explicitly written in (3.27) and (3.28) of [46]. 

While Cj and Cg receive the loop-level RS contributions calculated in Section 3, Cgy 
and C|q^ are corrected at tree level from the new flavor-changing couplings to the Z boson 
and the exchange of neutral electroweak gauge boson KK modes. In this analysis we only 
keep the leading contribution to each of these operators, i.e. we consider AC^^ g^ at one 

loop and ACgy^^^ at tree level. Strictly speaking, such an approach leads to an inconsistent 
perturbative expansion, but it is reasonable to expect that the one loop corrections to the 
latter Wilson coefficients are sub-dominant with respect to the tree level contributions, and 
by considering only the RS tree level contribution one should still capture the dominant 
NP effects. 

Explicit expressions for ACgy and ACJq^ can be straightforwardly obtained from [22]. 
These expressions can be written in terms of RG invariants AY^'' and AZ^'' and the 
coupling a, which itself is only very weakly scale dependent above Mw- Thus one may use 
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these expressions to directly write the RS contributions at the scale Mw, 

ACgv 



a 
2^ 



AY, 

' 4AZ., 



a 

2^ 



sin^ 6w 
AY' 



4AZ' 



sin^ 9w 
a An 
27r sin^ 9w 



(5.5) 
(5.6) 



^^lOA = -iTZ—2h (5-7) 



a 



AY' 



ACloA = -1^-^^ (5.8) 



27r sin^ 9w 
The functions Ay''-* and AZ''^ are given by 

-;^-^E^^%f^A^^(-), (-) 

^^^= v4^E ,,/F^^L A^-(X). (5.12) 

^ts^tb ^ 8M^g^y^ sm^ 9w 

Here the summation runs over X = Z, Z^^' ,A^^> in the minimal model and over X = 
Z, Zh, Z', A^^> in the custodial model. The flavor violating 4D fermion gauge boson cou- 
plings A^^(X) depend on the overlap of the fermion profile with the corresponding gauge 
boson profile. Their explicit form depends on both the fermion and gauge boson mixing 
matrices. The explicit expressions are complicated and unilluminating, hence we do not 
quote them here but refer the reader to appendix A of [22]. Furthermore 

V 2 27r sm 9w 

The tree level contributions to 6 — )• sfj,~^fi~ transitions in the minimal RS model are 
evaluated in [28] without making the approximations of taking into account only the first 
KK modes or treating the Higgs vacuum expectation value as a perturbation. In this paper 
we are mainly interested in the effects of ~ 2.5 TeV KK modes. As these are ruled out in 
the minimal model by precision electroweak constraints, we focus on the phenomenological 
effects of the custodial RS model on these transitions. 

For the study of observables related to 6 — t- sfi^fi^ , it is useful to introduce the effective 
Wilson coefficients at the scale fif, that include the effects of operator mixing, 

Cf = (Cf )sM + ACrifXb) , C'f = {C'f)sM + AC'^ifXb) , (5.14) 

C^W) = iCI^)sM{q') + -AC,v , C',f = -AC^y , (5.15) 

a a 

2tt 27r 

CtoA = (ClOyl)sM H ACioA , CfoA = — ^C'lOA • (5-16) 

a a 
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The SM values of the effective Wilson coefficients can be found in Table 2 of [64], which 
also gives the c^ dependence of (C|y )sm('?^) in terms of a linear combination of the other 
Wilson coefficients. While in principle all contributions have to be taken at the scale /i^, 
the NP contributions to C'gyioA ^^^ invariant under renormalization group evolution. 

With these effective Wilson coefficients at the B scale, we are now equipped to study 
observables in 6 — )• sfj,~^fi~ transitions. While this system offers a plethora of observables 
for study, a detailed analysis of all of them is beyond the scope of this paper, and we 
concentrate on studying a few benchmark observables that are particularly relevant for RS 
physics. A numerical analysis is presented in Section 6. 

In passing we would like to remark on the pattern of contributions to C'gyioA ^^^ ^^^ 
custodial model, as pointed out in [22]. Due to the suppression of flavor violating ZcTj^cfj^ 
couplings by the discrete Plr symmetry, the main contributions arise in the primed Wilson 
coefficients Cgy^^^, which are absent in the SM. Since the right-handed b quark, localized 
significantly further away from the IR brane than the left-handed one, is far less sensitive to 
flavor violating effects introduced by the RS KK modes, the RS effects in Y} , Zg turn out 
to be rather small (typically below 10%). This pattern is very different from the minimal 
model, where the Plr suppression mechanism is absent and large tree level flavor violating 
Z couplings to left-handed down-type quarks are present. 

5.2 Benchmark observables 
5.2.1 Br{B -^ XsfJ.+fJ.-) 

For very low lepton invariant mass q^ — )• 0, the B — )• Xgfi^n' transition is completely 
dominated by the photon pole and doesn't provide any new insight with respect to the 
B —7- Xs^ decay discussed above. Furthermore, in the intermediate region 6 GeV < q^ < 
14.4 GeV the sensitivity to NP is very small, as the decay rate in this region is completely 
dominated by charm resonances. Hence one usually restricts oneself to either the low q^ 
region 1 GeV^ < q^ < 6GeV^, or the high q"^ region g^ > 14.4 GeV^. In what follows we 
will consider only the low q^ region. While the high q'^ region is potentially interesting 
since it exhibits a small tension between SM prediction [65] and experimental data [66, 67], 
it is far less sensitive to NP in C^ , which is the main focus of this study. In the custodial 
RS model, the tension in the high q'^ region cannot be resolved since the new contributions 
to Cgy^Qy^ are generally small [22]. In addition, the high q'^ region is subject to larger 
theoretical uncertainties. 

In the low q'^ region, adapting the formulae of [68] to the more general case of complex 
NP contributions, we find 

Br(B ^ X,/i+/x-)'°" «' = Bt{B -^ Xsfi+fi-)sM '' + ABr(5 ^ X,;U+/i-)i°" "' (5.17) 

with the NNLL prediction [69] 

Br(B -^ Xsfi+fi~)sM ^' = (15.9 ± 1.1) • 10"^ (5.18) 
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and the NP contribution [68] 



ABr{B -^ XsH^fi 



+ //^^lo™ 9^ 



~ 10 



-7 



- 0.517 Re(AC7(/Ufc)) - 0.680 Re(AC^(/Xb)) 

+ 2.663 Re(<5C9y) - 4.679 Re((5CioA) 

+ 27.776 (|AC7(A.5)|' + |AC^(/ife)|2) 

+ 0.534 (I^Cgyp + IJC^v-P) 
+ OM3{\6Cioa\^ + \6C[oa\'') 

+ 4.920 Re {AC7inb)SC^v + AC;(;Ub)<5C7^V) 1 > (5-19) 



where we defined 



6Q. 



2n 



a 



ACi^ 



(5.20) 



Note that we dropped all interference terms between unprimed and primed contributions 
since they are suppressed by a factor ms/mb and therefore small. The only exception is the 
term linear in ACy, which receives a large numerical enhancement factor, and is therefore 
non-negligible; hence we keep it in our analysis. 

The measurements of BaBar [66] and Belle [67] yield the averaged value 



Br{B -^ X,^+/i-)|,°^'?' = (16.3 ± 5.0) • 10" 



(5.21) 



As LHCb is not well suited for performing inclusive measurements, a significant reduction 
of uncertainties will only be feasible at the next generation B factories Belle-11 and SuperB 

[70-73]. 

5.2.2 B -^ K°*(^ 7rK)/x+/x- 

While the inclusive B — )• Xsfi^fi'^ mode is theoretically very clean, such measurements are 
experimentally challenging, and competitive results (in particular for angular distributions) 
will not be available before the Belle II and SuperB era [70-73]. For this reason, exclusive 
decay modes have received well-deserved attention. An especially interesting decay is 
B — )• K*{-^ KTT)fj,~^n~ , where a plethora of angular observables can be studied thanks to 
the four-body final state [61, 64, 74-79]. These can provide detailed information on the 
operator and flavor structure of the underlying NP scenario. 

The downside is that many B — )• i^*(— )■ Ktt)^^ jjl^ observables, such as the branching 
ratio and differential decay distribution, are plagued by large theoretical uncertainties in the 
determination of the B — )■ K* matrix elements governed by long-distance non-perturbative 
QCD dynamics. These matrix elements are most conveniently described by a set of seven 
form factors. Presently, the best predictions for these form factors at large final state meson 
K* energies, i.e. small lepton invariant mass (f , stem from QCD sum rules at the light cone 
[80]. Furthermore, non-factorizable corrections are calculated using QCD factorization, 
which is only valid in the low q^ regime.^ On the other hand, as mentioned above, at very 

^Significant progress has recently been made on the form factor predictions in the large q^ region [81-84]; 
nevertheless we will not consider this kinematic regime since it is less sensitive to NP entering C\ than 
the low q^ region. 
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low q^ < 1 GeV^ the b — )• sfi'^fj.~ transition is dominated by the C^ contributions due to 
the infrared photon pole and therefore does not provide any insight beyond what is already 
obtained from b — )• 57. Consequently, we henceforth restrict our attention to the range 
IGeV^ <g2 <6GeV2. 

Fortunately, it is possible to partly circumvent the theoretical uncertainties by studying 
angular observables that are less dependent on the form factors in question. Detailed 
analyses of their NP sensitivity and discovery potential have been performed by various 
groups, both model-independently and within specific NP scenarios [61, 64, 74-77]. We 
leave such a detailed analysis in the context of RS models for future work. We focus 

instead on two benchmark observables, the forward backward asymmetry Ap^, which is 

(2) 
experimentally well constrained, and the transverse asymmetry Aj,, which offers unique 

sensitivity to NP in the primed Wilson coefficients. 

We note that the recently measured CP asymmetry Ag [85, 86], as defined in [64, 87], is 

also very sensitive to NP in Cj and therefore is in principle an interesting observable to look 

for RS effects. Because it is sensitive to the phase of CL it yields partly complementary 

information with respect to the CP conserving transverse asymmetry A)p . Although this 

CP asymmetry is theoretically very clean, contrary to those studied in [76], we leave a 

detailed study within RS for future work. 

Forward backward asymmetry The forward-backward asymmetry ApB in i? — )■ K*fi^iJ,' 
decays is defined by 

where 9^ is the angle between the K* momentum and the relative momentum of /i^ and 
/i~. ^FB has recently received a lot of attention as data from BaBar, Belle, and the 
Tevatron seem to indicate a deviation from the SM, albeit with low statistical significance 
[85, 88, 89]. On the other hand, recent LHCb data [90] show excellent agreement with 
the SM prediction, and as uncertainties are presently dominated by statistics, an improved 
measurement should be available soon. 

A precise theoretical determination of ^fb is appealing since it offers a sensitive probe 
of the helicity of NP contributions. To leading order, the forward backward asymmetry is 
proportional to [87] 

t2 \ / o»^2 



A-p-Q{q ) oc Re 



\Cw{q^) + "^'^7 j C'lOA - [C'w + "o2^^7 j c\ 



7* 



(5.23) 



where we dropped the superscript "eff" for the effective Wilson coefficients at the scale ^f,, 
(5.14-5.16). From (5.23) we can see explicitly that ^fb does not receive contributions from 
the interference of different chirality operators (unprimed and primed) . Consequently, with 
the SM contribution being the dominant effect, potential non-standard effects in A-pB arise 
mainly from NP in C7 and Cgy. On the other hand, ^fb is rather insensitive to NP in the 
primed Wilson coefficients C*/-,, gyioA- 

^FB has been studied in the context of the minimal RS model considering only tree 
level contributions and omitting loop level dipole contributions to C^ [28] , where small 
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positive contributions to ^fb were found. While ApB is very sensitive to NP effects in C7, 
the RS dipole contributions we calculated predict rather small contributions to this Wilson 
coefficient. On the other hand, ^fb is insensitive to Cy, where RS effects are expected 
to be more pronounced over the SM. Thus the overall prediction of small deviations of 
ApB from the SM obtained in [28] remains consistent with our calculations. Note that the 
restriction to tree level RS effects is not necessarily a good approximation for observables 
sensitive to C^, such as F^, which was also studied in [28]. A detailed study including 
one-loop contributions to the dipole operators would therefore be desirable but lies beyond 
the scope of the present analysis. 

In the custodial RS model, due to the protection of the Zd]^d-'^ vertex [25], the RS 
contributions to Cgv,ioA are highly suppressed, and only the new contributions to the 
primed operators are relevant. As ^fb is insensitive to the latter Wilson coefficients, it 
remains very close to the SM prediction. 

We conclude that RS effects in the forward backward asymmetry Afb are generally 
small, so the recent data from LHCb do not pose any stringent constraint on the minimal 
or custodial model, the latter being even more insensitive to RS contributions. 

Transverse asymmetry A^ The asymmetries Aj^, which are introduced in [75, 91], 
offer a particularly good probe of NP in 6 — t- s^i^ ^^ transitions since at leading order they 
are free of any hadronic uncertainties and are given in terms of calculable short distance 
physics. In this paper we will restrict ourselves to the study of the asymmetry 

A)p = ^-^^ ^ . (5.24) 

Here A±_ and A» are the transversity amplitudes [91] describing the polarization of the K* 
and the /^"^/x" pair; both are transverse with linear polarization vectors perpendicular (_L) 
or parallel (||) to each other. In the limit of heavy quark (m^ — )• cxd) mass and large K* 
energy (small q^), this asymmetry takes a particularly simple form [76] 

^ ^^ ^ \CioA? + \C[^a\'' + F^ (IC7P + |C;|2) + ICgvP + 2FRe(C7C9V) ^ ' ^ 

with F = 2mbmB/q'^, and we have again dropped the superscript "eff" from the Wilson 

(2) 
coefficients. In this limit it is clear that Aj. is independent of form factors and is governed 

only by calculable short distance physics, making this observable theoretically clean. Sec- 
ond, we notice that since the primed Wilson coefficients are highly suppressed in the SM, 

(2^ (2^ 

{A}p')gy[ is very small. Aj^' therefore offers unique sensitivity to NP entering dominantly 
in the primed operators C'r^-ygyioA- This asymmetry is thus a benchmark observable for 
discovering RS physics in B — )• K* fi^ fi^ decays. We investigate the possible size of RS 
contributions to this channel in our numerical analysis in the next section. 

A first measurement of Aj, by CDF [85] is still plagued by large uncertainties. LHCb 
has recently put more stringent constraints on this asymmetry, and more precise measure- 
ments will be possible in the near future [86]. With 10fb~ of data, LHCb is expected to 
reach a sensitivity of about ±0.16. 
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6 Numerical analysis 

6.1 Strategy 

In this section we present a numerical analysis of the observables introduced in the previous 
sections. To this end we fohow the fohowing strategy: 

1. The first goal is to understand the generic pattern of effects induced by RS penguins 
on flavor observables. We generate a set of parameter points that satisfy the known 
experimental constraints from quark masses and CKM parameters. However, we do 
not yet impose any additional flavor bounds so as not to be biased by their impact. 
With these points we evaluate the new RS contributions to the Wilson coefficients 
ACy and ACg at the KK scale for both the minimal and the custodial model. 
Subsequently we calculate the new contributions to the branching ratios of B -^ ^s,dl 
and analyze the constraints. 

Note that the same set of parameter points is used for the minimal and the custodial 
model in this case, in order to minimize the sampling bias on the results obtained. 

2. The second goal is to understand the effect of the RS penguins on the existing pa- 
rameter space for realistic RS models. We restrict our attention to the custodial 
model, which can be made consistent with electroweak precision tests for KK scales 
as low as Mkk — 2.5 TeV. In addition to quark masses and CKM parameters, we 
now also impose constraints from AF = 2 observables which are analyzed at length 
in [25] . After evaluating the size of the effects in the B — )• Xg^dl branching ratios and 
their constraint on the model, we study the benchmark observables outlined above, 
namely the CP asymmetry in S — )■ K*^^ the branching ratio Br(i? — >• XsH~^fi~), and 

,(2) 



frf\ 

the transverse asymmetry Aj, in S — t- K*fi^fi'' decays. 



Throughout our analysis we restrict ourselves to 1/R' = 1 TeV, so that the lowest KK gauge 
bosons have a mass of Mkk — 2.5 TeV. We note that in the minimal model such low KK 
masses are already excluded due to unacceptably large corrections to electroweak precision 
observables. However, we use the same mass scale for both the minimal and custodial 
models to enable a straightforward comparison of the two sets of results. Furthermore, we 
restrict the fundamental Yukawa couplings to lie in their perturbative regime, i. e. \Yij\ < 3. 
More details on the parameter scan can be found in [25]. 

6.2 General pattern of RS contributions 

This part of the numerical analysis is dedicated to determining the size of NP effects 
generated by the RS KK modes in the dipole operators Cj, Cj and Cg, Cg mediating the 
b — )• (s, d)'y and b — )• (s, d)g transitions respectively. We advise caution when interpreting 
the density of points since these distributions are influenced by the details of the parameter 
scan performed. The qualitative features in our plots should however remain unaffected by 
the scanning procedures. 
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Figure 4. RS contributions to the fe — > s Wilson coefficients Cr{MKK) (upper left), C''^{Mkk) 
(upper right), Cs{Mkk) (lower left) and Cg(MKK) (lower right) in the minimal (red, dashed) and 
custodial (blue, solid) models, and from the misalignment contribution alone (black, dotted). 

The first row of Fig. 4 shows the RS contributions to C7(Mkk) and Cj^M^k) in the 
6 — )• s system. Observe that the total RS contribution (dashed red and solid blue distribu- 
tions, corresponding to the minimal and custodial model) to the primed Wilson coefficient 
is typically an order of magnitude larger than the corresponding effect in the unprimed 
Wilson coefficient. This matches the naive expectation that the bi — )• sr transition should 
be enhanced relative to 6r — >• sl due to the hierarchy /q^ » /^^ of fermion localizations. 
Furthermore the custodial contribution is somewhat enhanced relative to the minimal one, 
due to the additional fermion modes running in the loop. Also shown, in black (dotted), 
is the contribution to C7(Mkk) and (^^(Mkk) generated by only the misalignment term, 
which is equal for the minimal and the custodial models. Unlike the anarchic term, this 
contribution is generically comparable in both cases. This naively unexpected behavior 
is explained in Appendix C. While it is subdominant but non-negligible in the case of 
C7(Mkk)) it turns out to be generally irrelevant in the case of Cj{Mkk)- 

The second row of Fig. 4 shows the results for the gluonic penguin Wilson coefficients 
Cs and Cg. The values at the KK scale are larger than the corresponding values of Cj 
and C'y by about an order of magnitude due to the large contribution from the diagram 
containing the non-Abelian SU{3)c vertex, which is absent in the b ^ s'j penguin. Other 
than that, the pattern of effects is qualitatively similar to that for Ci : the primed Wilson 
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Figure 5. RS contributions to the b —>■ sj Wilson coefficients Cy (left) and Cy (right), evaluated at 
the scale ^{, = 2.5 GeV. The minimal model distribution is shown in red (dashed), and the custodial 
one in blue (solid). 
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Figure 6. Relative sizes of anarchic contributions to the Wilson coefficients C7(/i(,) (left) and 
C'i{^b) (right) from the RG evolution and operator mixing of ACg from A/kk to /ifc, normalized to 
the Higgs penguin contribution to ACy (Mkk)j with relevant RG evolution factors included. The 
black (dotted) peak shows the ratio of the Higgs penguin contribution to ACg (Mkk)- The red 
(dashed) and blue (sohd) distributions show the ratio of the gluon penguin to ACg (Mkk) for the 
minimal and custodial model respectively. 



coefficient is larger than the unprimed coefficient by about an order of magnitude, and 
the custodial model yields somewhat bigger effects than the minimal model. Furthermore, 
the misalignment contributions to the unprimed and primed Wilson coefficients are again 
roughly comparable; consequently, its effect is negligible in Cg but can be sizable in C%. 

To facilitate comparison with other models of NP, Fig. 5 shows the RS contributions 
to the 6 — )• S7 Wilson coefficients C7 (left) and Cy (right) evaluated at the scale /i^ = 
2.5 GeV, i.e. taking into account the RG evolution and operator mixing with Cg . The RS 
contribution to Cy turns out to be small and typically constitues less than a few percent of 
the SM value Cy(//fc)^'^ = —0.353. On the other hand, Cy is suppressed by ms/mi, in the 
SM, so the unsuppressed contribution from RS dominates, though its value is still typically 
smaller than Cy(/ib)^^. 

Next, we examine the relative importance of the various RS contributions to the ef- 
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fective b ^ s'j Wilson coefficients at the scale //;,. Fig. 6 shows the size of the two main 
anarchic contributions to ACg (Mkk) (see Fig. 2a for the relevant Feynman diagrams) nor- 
malized to the anarchic contribution to AC^ (Mkk) (see Fig. la). For a straightforward 
comparison, we also include the relevant RG evolution factors from eq. (4.6). The ratio 
of the Higgs penguin contribution to AC^ (Mkk) and ACg (Mkk)) shown by the black 
(dotted) peak, is constant and equal for both the minimal and custodial model. As the 
relevant diagrams depend on the same loop integral and the same combination of Yukawa 
couplings, their relative size at the KK scale is simply given by the electric charge Qu of 
the up-type quark coupled to the photon. After including the RG running down to the 
scale fib, the Higgs penguin contribution to Cg turns out to be roughly a 50% correction 
to the effect of the anarchic ACy (Mkk) contribution. 

The effect of the gluon penguin diagram in ACg (Mkk) depends on a different loop 
integral and a different combination of Yukawa couplings than the Higgs diagram in 
AC7 (Mkk)- Consequently its relative size, again including the relevant RG factors, varies 
considerably within the minimal (shown in red, dashed line) and the custodial (shown in 
blue, solid line) model. Observe that the distribution for the minimal model is rather 
symmetric and peaked around 1, implying that the RS b ^ sg loop generally contributes 
as much as the RS 6—7-57 loop in low energy observables, even yielding the dominant RS 
contribution in parts of the parameter space. This is in contrast with the SM case, where 
the Cg contribution only gives a few percent correction to the dominant Cj contribution. 
In the custodial model the gluon penguin contribution becomes even more important, so 
that the peak of the distribution gets shifted above 1. Since, as opposed to the Higgs 
penguin, the additional custodial gluon penguin diagram shown in Fig. 3 carries the same 
Yukawa spurion as the minimal model diagram, they simply add constructively, further 
enhancing the effect of the gluonic penguin contribution. Neglecting these contributions or 
even the Cg contribution as a whole, as sometimes done in the literature, would therefore 
be a rather poor approximation. Note that the relative importance of the gluon penguin 
diagrams depends crucially on the matching of the 5D to the 4D strong gauge coupling. 
Invoking one loop level matching rather than tree level matching as done here whould 
reduce their relative size by roughly a factor of four. On the other hand the presence of 
brane kinetic terms could further enhance the gluonic penguin contribution. 

For the sake of completeness Fig. 7 shows the Wilson coefficients for the 6 — )- d system, 
in analogy to Fig. 4. The pattern of effects is very similar to the case of the 6 — )• s system 
discussed above. 

Fig. 8 shows the predicted deviations from the SM in the B — )- X^.d^ branching ratios 
in the minimal and custodial models. We observe that in both models these branching 
ratios typically obtain a moderate positive NP contribution well within the current exper- 
imental and theoretical uncertainties. Nevertheless, the decays in question put nontrivial 
constraints on parts of the RS parameter space and should be included in a complete anal- 
ysis of RS flavor phenomenology. As expected from the size of the Wilson coefficients, the 
custodial model induces somewhat larger effects than the minimal model. 

Interestingly, this pattern of effects is very different from that of the ADD model 
of a universal extra dimension [92], where the KK excitations affect mainly the Wilson 
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Figure 7. RS contributions to the h ^ d Wilson coefficients C7(AfKK) (upper left), C7(A/kk) 
(upper right), C8(AfKK) (lower left) and Cg(AfKK) (lower right) in the minimal model (red, dashed), 
the custodial model (blue, solid), and ftom the misalignment contribution alone (black, dotted). 

coefficient Cy, while the opposite-chirality Wilson coefficient C^ remains very small [93, 94]. 
Since the ADD contribution interferes destructively with the SM contribution, a rather 
pronounced suppression of Br(i? — )• ^^7) is predicted, which was used in [95] to derive the 
bound 1/R > 600 GeV on the radius R of the extra dimension. 

We also investigated the dependence of the size of the RS contribution to Br(i? — )• Xsj) 
and {Bi{B — )• Xd'y)) on the average Yukawa coupling y^,, but did not find any significant 
correlation within our parameter scan. These findings at ffist sight seem to contradict the 
analytic estimate in section 4.4. Recall that these estimates have been performed in the 
fully anarchic limit where Kk is the only free parameter in the Yukawa sector. On the 
other hand, our scan varies all independent parameters in the flavor sector, so that 0(1) 
deviations from the fully anarchic ansatz are intrinsic. The dependence on these additional 
parameters fully hides the dependence on y*; note also that the latter only varies over an 
0(1) range. 

6.3 Effects on benchmark observables 

We now restrict our attention to the custodial model and consider only parameter points 
that agree with the existing constraints from AF = 2 transitions, as analyzed in [25]. We 
also impose the bounds from the B — )■ Xg^al decays as approximated in (4.3-4.4), so that 
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Figure 8. RS contribution to Bi{B -> Xg'^) (left) and (Br(i? — ?> X^'y)) in the minimal (red, 
dashed) and custodial (blue, solid) model. The experimental constraints according to (4.3) and 
(4.4) are displayed as grey bands. 
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Figure 9. CP asymmetry Sk'-^ as a function of Br(i3 — >■ ^^7). The black-and-white dot indicates 
the central SM prediction, while the dashed lines show the experimental central values. The grey 
bands display the experimental Ict and 2(7 ranges for Sk*^- We omit showing the uncertainty in 
Br(i3 -^ Xg^) as it covers the whole range. 

all points displayed in the plots lie within the experimentally allowed region. 

Since the dipole operators depend on a different combination of RS flavor parameters 
from the tree level contributions to AF = 2 processes [25] and /S.F = 1 rare decays [22], 
observables related to the various sectors are essentially uncorrelated; hence we do not 
show any numerical results here. 

Fig. 9 shows the correlation between the time-dependent CP asymmetry S'x*7 and the 
branching ratio of B — t- Xg'y. Observe that Sk'-^ can receive large enhancements relative 
to its tiny SM value. While non-standard effects in 5x*7 are possible for any value of 
Br(S — )• Xs7), large effects are more likely with enhanced values of the branching ratio. 
This is related to the fact that RS contributions dominantly affect C'-j. While the SM 
prediction for B — t- Xg^ is in good agreement with data, it lies below the central value. 



26 



r 1.60 



•.%^. •«• 







.*«'••• 




3.2 3.4 

10'' 61(3^X53-) 



Figure 10. Correlation between Br{B -> Xg-f) and Br(i3 -> XsfJ.^n~) for g^ e [1,6] GeV^. 
The black-and-white dot indicates the central SM prediction, while the dashed lines show the 
experimental central values. We omit showing the experimental and theoretical uncertainties as 
they cover the whole range. 




Figure 11. Transverse asymmetry Arp' as a function of q^, for a few parameter points. The SM 
prediction is indicated by the thick black line, while each blue line corresponds to an RS parameter 
point. 

and an enhancement of this branching ratio is preferred. One can also see that large 
enhancements are possible in SK*'y, and that the present experimental 2a range excludes 
only a small fraction of the RS parameter space. 

The decay B — t- Xsfi~^^~ poses strong constraints on various extensions of the SM, 
hence it is worth studying it in the custodial RS model. Fig. 10 shows the custodial RS 
branching ratio Br(i? — )• Xg^i^ ^i~) in the low q^ region as a function of Br(i? — )• Xg'^)- We 
observe that the enhancement in the custodial RS model is rather small, typically below 
10%. Due to the experimental and theoretical uncertainties involved, this channel does not 
put any significant constraint on the model. 

Observables far more sensitive to NP in C'-j can be constructed from the angular distri- 
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Figure 12. Correlation between Sfc'-y and A^ 



i(^'(g2 ^ IGeV^). The black-and-white dot indicates 
the central SM prediction, while the dashed line shows the experimental central value. The grey 
bands display the experimental la and 2a ranges for Sk'j- 



bution of i? — )• K*fi~^fi~. Of particular interest is the transverse asymmetry ^4^, , whose q"^ 
dependence is shown in Fig. 11. Observe that large enhancements relative to the small SM 
value are possible, in particular in the very small q^ region < 2 GeV . This pattern can be 
understood from (5.25): the (74 contribution is enhanced at small q'^ due to a 1/g^ factor, 

(2) 

see also [76, 96]. Also note that the custodial RS model predicts a zero crossing for A}p 
at q^ ~ 2.7 GeV . The differential asymmetry would exhibit a very different shape if the 
dominant NP contribution appeared in C^qa- This underlines the model-discriminating 

(2) 

power of the A)p asymmetry — in the custodial RS model a deviation from the SM is most 
likely to be observed for small q'^, whereas other models that dominantly affect C(q^ predict 
large effects for larger q^. This pattern is particularly interesting in light of LHCb and the 
next generation B factories, which will soon be able to measure this asymmetry. 

(2) 

Finally, one may consider a possible correlation between Sx'-y and A}p . Both observ- 
ables are mostly affected by a large CL hence some nontrivial correlation can be expected. 

(2) 

On the other hand, SK*-y is CP violating while Aj^ is CP conserving, so the phase of Cy 
can wash out such correlations. Fig. 12 shows A)p {q"^ = IGeV^) as a function of Sx'-y, 
where a nontrivial linear anti-correlation is seen between the two observables in question. 
However, this correlation is visibly weakened by the impact of the phase of C^, as expected. 



7 Conclusions 

In this paper we have performed an explicit 5D calculation of the dominant contributions 
to the Wilson coefficients C7, C7, and Cg, Cg that mediate the b ^- Sjd'j and b ^ s,dg 
transitions respectively, in the RS setup with bulk fermions and gauge bosons and an 
IR-brane localized Higgs. We have evaluated the relevant diagrams for both the minimal 
scenario with only the SM gauge group in the bulk, and for the custodial model with the 
electroweak gauge group extended by SU(2)r and a discrete Plr symmetry. Our main 
findings from this analysis can be summarized as follows: 
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• The RS contributions to Cj typically exceed those to Cy by an order of magnitude, 
and the latter remain a rather small correction to the SM value. This pattern can be 
understood by considering the bulk profiles of the quark fields involved: the primed 
Wilson coefficient describes the decay of a left-handed h quark, which, being localized 
towards the IR brane, is more sensitive to flavor violating effects than the right- 
handed h quark entering C7. Analogous comments apply regarding the hierarchy 

• Contrary to the SM, where Cs < C7, RS contributions to the gluonic penguins are 
larger than the ones to the photonic penguins. This results from the large contri- 
butions from the diagram containing the non-abelian triple gluon (KK gluon) ver- 
tex, which is absent in C\ and does not change flavor in the SM. In addition, the 
renormalization group mixing of Cy and Cg is more pronounced due to the large 
separation of the Mkk and mf, scales. Consequently, gluonic penguin contributions 
have a significant impact on h — )• 5,^7, comparable to or larger than the photonic 
penguin contribution. This is in contrast to the SM, where they yield only a few 
percent correction to the photonic Wilson coefficients at the rrn, scale. 

• In all cases, the dominant effect comes from the anarchic contributions, which are 
not aligned with the SM quark mass matrices. However, the unprimed (right to 
left) operators pick up appreciable contributions from misalignment diagrams, which 
are proportional to the SM quark mass matrices up to a dependence on the bulk 
spectrum. This is because, in contrast to the anarchic diagrams, the misalignment 
diagrams are not suppressed by the bn wave function relative to the bi wavefunction, 
as explained in Appendix C. 

• The impact on the Wilson coefficients in question is somewhat larger in the custodial 
model than in the minimal model, since the extended fermion content that was intro- 
duced to reconcile the model with the Zbb constraint yields additional contributions. 

For a study of the phenomenological implications of these new contributions, we re- 
stricted our attention to the custodial model since the minimal model is not consistent 
with electroweak precision constraints for low KK masses Mkk = 2.5 TeV. To this end, 
following [25] we performed a parameter scan of the 5D bulk masses and fundamental 
Yukawa coupling matrices, imposing constraints from quark masses and CKM parameters 
and from meson-antimeson mixing. We studied the bounds provided by the branching 
ratios Br(i? — )• Xg^) and (Br(i? — )• X^'j)) and the effects in a number of benchmark ob- 
servables, namely the time-dependent CP asymmetry Sk'-j, the inclusive branching ratio 
Br(i? — )• XsfJ-^fJ-~) and the forward-backward asymmetry ApB and the transverse asym- 
metry AIj, in i? — )• K*fi'^fi~, where we found the following patterns: 

• The branching ratios of the radiative inclusive B — t- Xg^^l decays provide a non- 
negligible constraint on RS models and exclude roughly 15% of the parameter points 
generated for the custodial model that were in agreement with bounds from AF = 2 
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observables. A complete phenomenological study should therefore take these con- 
straints into account. However, since the major part of parameter space survives, no 
useful bound on the KK scale can be derived. 

Due to more precise data and SM theory prediction, Bv{B — )■ Xg^) generally puts 
a stronger constraint on the RS parameter space than (Bt{B — )■ Xdj))- The lat- 
ter observable is still useful as it yields complementary information on the allowed 
parameter space. 

• As the RS contributions enter dominantly through the primed operators, a modest 
enhancement of the B — )■ Xg^d'y branching ratios can be expected, although a slight 
suppression is not rigorously excluded. Such an enhancement would be welcome in 
B — 7- Xs7, where the data lie somewhat above the SM value, albeit still in good 
agreement. On the other hand, for B — )■ X^'j the central values of the SM and the 
data are in excellent agreement and the uncertainties are sizable, and no prefered 
sign for the NP contribution can be deduced. 

• The inclusive branching ratio Br{B — )■ XgH^fi^) and the forward backward asymme- 
try ^FB in i? — )• K* jji^ iJL~ receive very small corrections from RS physics and remain 
in good agreement with recent data. While we restricted our analysis to the low q^ 
region, these statements also apply to the high q^ region since the latter region is 
mostly sensitive to NP in the electroweak Wilson coefficients C'gyioA' which remain 
SM-like in the custodial model. 

• We identify the time-dependent CP asymmetry S'x*7 in S — )■ K*^ decays and the 
transverse asymmetry Aj. in the low q^ region of S — t- K* ^i^ ^i^ as promising bench- 
mark observables to look for large effects generated by the custodial RS model. Both 
observables are known to be very sensitive to the primed Wilson coefficients, in par- 
ticular C7, which is dominantly affected by RS contributions. Furthermore, studying 
the (f' dependence of Aj, allows for a clear distinction of models such as the custodial 
RS model that dominantly affect Cj from models that predict large NP effects in the 
electroweak Wilson coefficient C^qa- 

In summary, our analysis shows that radiative and semileptonic B decays offer intrigu- 
ing possibilities to find deviations from the SM generated by RS KK modes and anarchic 
Yukawa structure. If such effects are found at the LHCb and the next generation B fac- 
tories, it will be particularly interesting to study the plethora of observables provided by 
these decay modes in a correlated manner, which offers the ability to distinguish RS with 
custodial symmetry from other NP scenarios that predict a different pattern of effects. 
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A Dimensionless Integrals for Leading Diagrams 

This appendix defines the dimensionless integrals associated with the leading contributions 
to the a and b terms of the dipole Wilson coefficients C^^s in Section 3. Details of the 
derivation of these integrals are found in the appendix of [39]. In the mass insertion 
approximation the Standard Model contribution appears as an infrared pole, which we 
subtract. 

A.l Propagator functions 

We use dimensionless integration variables x = kEZ G [wy, y] and y = kER' G [0, oo], where 
kE is the Euclidean loop momentum and w = (R/R') is the warp factor. The integrals 
are expressed with respect to the functions that appear in the mixed position-Euclidean 
momentum space fermion propagator, 

where the F functions are defined for x > x' (i.e. z > z') by 

~L ^ jxx'f/'^ Sc^ix-,y-)Sc^{x'_,wy-) -^ ^ {xx'f''^Tc^{x+,y-)T^^{x'^,wy-) 
y^ ScL{y-,wy-) + 2/5 Sci^[y-,wy-) 

(A.2) 

pR ^ jxx'f/'^ Tcj,{x^,y+)Tci,ix'^,wy+) ~^^ {xx'fl"^ Sc^{x+,y+)Sc^{x\,wy+) 
y^ ScR{y+,wy+) + y5 Scjiiy+,wy+) 

(A.3) 

The analogous functions for x < x' are given by replacing x -H- x' in the above formulas. 
S and T function are products of Bessel functions, 

Sc{x±,x'^) = 4±i/2(2;)i^c±i/2(a:^') - Ic±i/2ix')K^^i/2{x) (A.4) 

Sc{x±,x'^) = Ic±i/2{x)Kc^i/2{x') - Ic^i/2{x')Kc±i/2{x) (A. 5) 

Tcix±,x':^) = Ic±l/2{x)Kc^i/2{x') + Ic^i/2{x')Kc±i/2{x) . (A.6) 
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Similarly, the mixed position-Euclidean momentum space vector propagators are —i'q^'^G 
and iG for the 4- vector and scalar parts respectively. For x < x' , the G functions are, 



(^^ (x x') = (R')^^^' Tio{x,y)Tio{x', wy) 
R ^ ' R y Soo{wy,y) 

{R'?r^ ( M _ {RTxx' Soo{x,y)Soo{x', wy) 
R ^^^"^'^^- R y Soo{wy,y) 



G,{z,z') = '-^Gy{x,x') = '-^— '''T'. ' . ' ^ (A.7) 



where 



Tijix, y) = h{x)Kj{y) + Ij{y)K,{x) (A.9) 

Sij{x, y) = h{x)Kj{y) - Ij{y)K,{x). (A.IO) 

For z < z' the above formula is modified by x -H- x'. 

A. 2 Cr integrals 

We label vertices such that the external fermion legs attach to vertices 1 and 3, and the 
photon or gluon is emitted at vertex 2. Propagators attached to the brane x = y signify 
Yukawa couplings or mass insertions, which may change the fermion flavor as labeled by 
its bulk mass, c. We have left this c dependence implicit in the following expressions. 



Ic7a = dydxy 

\x 



2 fy\ [ _ r,pLyx rpLxy pRyy ^ 

4 1 2 

^Lyx pLxy ^Ryy V f Pi j^^y^\ p^^J/ J^RVV 2/ 

+ ^+,y ^+,y ^-,y 7^^^77Z^157^2^2 " o l^ ^^^^+,2/ J ^+,y ^- 



(y2 + {MwR'W 2 V ^^ +'y J +'y -^y y^ + {MwR'? 



2 V!'*^°-^-r) °-^i? ^"T y^ ^ iM„Rr + ^^« °-^-? °-^» %^ + (M>.i.p 



_l_ pLyy pRyx pRxy V \--fiif) F^yy\ p^y^ pR^y V 

^ ^+,y ^-,y ^-,y y2 + [MwR'Y 2 V^^^^+^y ) -^y -^y y2 + {MwR'y 



TpLyy ( n T7Ryy\ rpRxy V 



+ -FX%'[yd,,F-y^)F 



2 



'y y2 + [MwR'Y 



+ ^Fi,f (.*.«,?«-) J)-gy ^,^(,,^^^,), ] . (A.n) 

The C'jh integral is the sum of two parts corresponding to diagrams with an internal gluon 
(G) or scalar gluon (G5) in the loop, 

Each of these terms include diagrams with a single mass insertion, either on the incoming, 
internal, or outgoing fermion line. 
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/C76 = I dy dxi dx2 dxsy { —] dk^G^^ 

2+CL / ,, \ 4 



+ 2V^j Uj ^+^+,(-.i?') [D-P-,y P-,y + F+,y D-P-,y 

(\ 2+ci^ / \ 2—Cft 
^) (^) ( - ^+^f,f D-F^f p+f + y' ^^,f p^f p+S 



2 I X2 



2+cl / \ 4 

— ' r) jpLixsniiR' /y){mtR') 



^Cn = dy dxi dx2 dx^ - i 

(f^]^ D+F^f (y dk.Gf + 4 Gf) + y Gf{D+F^f d^J^]^ - F^f du.D^F^]^ ) 

.fvYfy V"* n pR(rn,R'){x,m,R'/y) 

(Ff 1/ 5+Ff/ (y a,,Gii +AGf) +yGf{D^F^f d^J^^ - F^f d^.D^F^]^ 



xi) V^s/ 



(Z?+F^,i,2 (4 + y 9.,)(F^,y i?+F^,f Gf) - y F^]^dk, (D+F^J D^F^f ) Gf 

+ D^F^^y D^F^' {4 +ydk,){F^fGf) - y D^F^ly F^^' Gf d^, D^F^f)} 

(A.14) 

A. 3 Cs integrals 

The Csa integral contains a piece identical to the C^a integral associated with the charged 
Higgs loop as well as gluon loop diagrams with three mass insertions, 

Ic3„=/ai + 24:i+/g„. (A.15) 

The gluon loops are labeled by the number of internal mass insertions, so that I^J is 
associated with the diagram with an external mass insertion on each leg, and the factor of 
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two on I^ accounts for the two possible placements of the external mass insertion^. 

C+F«»;D_i?«»D_ff»/D+Fi5{-|!,8,,(G;^Gf) + 10Oj^Gf}, (A.16) 

rpLly rpRyy j^ TpLy3 p. rpL{x3nn,R' /y){mt,R') r. .^12 ^23n, i K ^7\ 



2+CL / „, \ / „, \ 2-CH 



pLly pRy, ^Lyy ^Ryy | _ 5^ g^^ ^^12 ^23^) + 10^^,^^^}. (A.18) 



For Cgf,, the only dominant diagram is the gluon loop with an internal mass insertion. All 
other analogous diagrams (e.g. mass insertion on an external leg, or loops with G^) contain 
no zero modes and hence give negligible contributions after alignment. 



2+CL / ,, \ / „ \ 2-Cfl 



Ic,, = jdydx,dx,dx,{£^ (^)(^) 

'f { - \y dk, {Gf Gf) + 10 G"^ Gf]. (A.19) 
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B Charged Higgs diagram calculation 

As an example of how to calculate diagrams in the mixed position/momentum formalism, 
we present the calculation of the leading contribution to the anarchic piece of the Cy 
operator coming from the charged Higgs diagram in Fig. la. As discussed in Section 3.2, 
it is sufficient to compute the coefficient of the p^ term in the amplitude. This allows us 
to directly write the finite physical contribution to the amplitude without worrying about 
regularization of potentially divergent terms. In addition to the bulk fermion propagators 
in mixed position/momentum space, A(p,z,z'), which are given in Appendix A.l, the 
relevant Feynman rules are given by 

D Q f f 

H ^M 

A derivation of the propagators and a more complete set of Feynman rules is given the 
appendix of [39] . The amplitude for the diagram with a 6 of momentum p decaying into a 




These integrands differ by L -o- -R, but the integrals are approximately the same. 



34 



photon of momentum —q and a s of momentum p' is 



M^ 



J- f d'k 



R' /^x 4 



dz^^j UQ^{p')fQAG%,fD,nDM^H{k-p) 

(B.l) 



where k and £ index the flavors of the internal fermions and Ah is the 4D Higgs propagator. 
Writing k' = k + q, the Dirac structure G^ for the diagram with the mass insertion before 
(a) or after (b) the photon emisison is 



<-(a) 



ke 



ke 



AD,{k',R',z) Y AD,{k,z,R') AQ^{k,R',R'), 
An,{k',R',R') AQ^{k,R',z) 7^^ AQ^{k,z,R'). 



(B.2) 
(B.3) 



We may now expand the fermion propagators in terms of scalar functions F, which are the 
Minkowski space versions of the F functions defined in Appendix A.l to simplify the Dirac 
structure and write the integrand in the form 



^Qs (at^r^ + ati'^'r) Pr ud, Anik - p) 
where g^"^' is a scalar function that takes the form 



n e {a, 6} 



Sa) 



9ta^^k,k') 



F-(k',R',z) F-{k,z,R') F+(k,R',R' 



.W/ 



g'^>{z,k,k')= F^(k',R',R') D^FQ(k',R',z) D+F+{k,z,R' 



(B.4) 



(B.5) 
(B.6) 



The derivative operators D± are defined in (A.l). g^'"'' has a similar definition but, as we 
show below, drops out of the final expression. 

To identify thep^ coefficient, which in turn determines the coefficient of the C7 effective 
operator, Taylor expand in p and q and perform the integral. It is sufficient to take only the 
leading order terms since higher terms are suppressed by the ratio of the external fermion 
masses to the characteristic loop energy scale (e.g. niH or 1/R!). The terms proportional 
to g^"'' and g^^' thus can be expanded as 



g{z,k,k')$'^^'AH{k-p)= {g + 



k 



dk' k 



k'=k 



{it + <^)j'' {AHik) + 2k-pAJ,ik)) 



(B.7) 



The g terms yield expressions proportional to ■y^p and j^]^' ■ By using the Clifford algebra 
and the equations of motion for the external particles one can show that these terms are 
proportional to rrih'y^ and 2p''^ — rus'j^ respectively. Thus these terms can be ignored since 
these do not contribute to the p^ coefficient. The g terms, on the other hand, contribute 
expressions of the form 



I^PgAUk) - \k^^H{k) - 2gAH{k)^ p^ 



(B.^ 



where we write k = \fk^xh^ and g is evaluated at (; = 0, i.e. k' = k. 
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Finally, the coefficient Uki of the amplitude (3.5) can be written with respect to the 
Wick- rotated integral of the prefactor multiplying p^^, 

aM = -2^ Y: l^dyl\xy'[iyhg(-)Aj, + l^ + 2g(AAH, (B.9) 

where we have defined the dimensionless integration variables x = —ikz, y = —ikR', 
and y' = —i{k + q)R' ■ This is equivalent to replacing the Minkowski space functions 
F{kR' ,zR' ,z'R') with the Euclidean space functions F{y,x,x') defined in Appendix A.l. 
The g and Ah functions are evaluated at k ^ iy and rriH — )• mj{R'- These are now 
completely scalar expressions that can be evaluated numerically. The explicit form of the 
integrand is given in (A. 11). 

Other diagrams are calculated following a similar algorithm with the caveat that dia- 
grams with bulk gauge bosons have 5D propagators, which carry additional space integrals 
over the extra dimension. 

C Estimating the size of the misalignment contribution 

In this appendix we clarify a subtlety in the size of the anarchic contributions {ACj g^) 

versus the misalignment contributions (AC^ g^,) to the Wilson coefficients, as defined in 
Section 3.2. For the anarchic contributions the relative sizes of the right-to-left (unprimed) 
coefficients to the left-to-right (primed) coefficients are given by the relative size of the fb^ 
and fbj^ wavefunctions on the IR brane. On the other hand, the misalignment contributions 
for the two chiral transitions do not follow this pattern and are, in fact, of the same 
order of magnitude. We show here that this apparent inconsistency can be understood by 
accounting for cancelations coming from the rotation to the SM fermion mass basis. 

For simplicity, consider the 2x2 matrix of misalignment diagrams qf" — )• qf" where we 
only consider the second and third generations. This transition is given by the bij term in 
(3.5), which we may parameterize as 

/ . ,. ,, ^ I{b-c-d)yu{b-c + d)yi2\ .„^. 

(misalignment term) • ~ ., , ,. ., , .. . (C.l) 

■> y{b + c-d)y2i {b + c + d)y22j 

Here we have written b as an average scale for the bij matrix, and yij = fQiYj^^fD ■ The 
c '^ 10^^ and d ~ W~^ terms represent deviations from the average. In particular, the c 
deviations account for the effect of an internal bi (whose zero mode profile is very different 
from that of the light quarks) while the d deviations account for the smaller effect of an 
internal 6r. 

In order to pass to the physical basis, one must apply to this matrix the same rotation 
that diagonalizes the SM mass matrix, which is proportional to y. The off-diagonal terms 
of the rotated misalignment matrix give the Cj and Cj coefficients (the argument for Cg 
is identical), 

(misalignment term) [ , ,, ) ~ ( ^/f ) • (C2) 
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The parameters 5 and 7 are ratios of the left- and right-handed zero mode wavefunctions 
on the brane; the primed and unprimed parameters are related by a minus sign. 

We focus on order of magnitude estimates, so we introduce a numerical parameter 
e ^ 10^^. Normalizing the Yukawa to 2/22 = Ij our parameters are approximately 

c ~ e d ~ e yii ~ e yu ~ e 1/21 ~ e ^ ~ e 7 ~ £• (C3) 

Note that e is merely a fiducial quantity, not an expansion parameter of the model. We 
now apply the rotation (C.2) and study the order of magnitude of the off-diagonal terms. 
By construction the terms proportional to b are completely diagonalized. We consider the 
terms proportional to c (/fe^) and d (fbji) separately. 

C.l Misalignment from /b^ 



First consider the terms proportional to c, which are split by the relative size of fb^ versus 

JSL 

like 



fs^ from internal zero mode propagators. The part of the Cy|^ term proportional to c goes 



C: 



[y2i + 1 y22) - S' {1' yi2 + yn) ■ (C.4) 



Naively the first term is of 0{e) and appears to dominate the expression. This, however, 
does not account for relations coming from alignment. Observe that the minus sign here 
comes from the choice of parameterization in (C.l). Further, observe that changing the 
relative sign in (C.4) is equivalent to changing the sign of c in the top row of (C.l). In 
this case, however, the c matrix would be completely aligned with the SM mass matrix 
and the off diagonal term (C.4) would vanish. Thus the first and second terms in (C.4) 
must be of the same order of magnitude in order for them to cancel when the relative 
sign is swapped — in other words, (y2i + 7' 2/22) ~ e^ in order to match the naive order of 
magnitude of the second term. We thus have 

cC'jb\,^e^ (C.5) 

This observation reflects the key cancelation that causes the relative size of the primed and 
unprimed misalignment terms to differ from that of the anarchic terms of the amplitude. 
The contribution to the C-^b term proportional to c is 

Cjblc'^ Hi 2/21 + ^22) - (7 yii + yi2) • (C.6) 

Unlike C^^^, both terms in the above expression are dominated by their 0{e'^) components 
and we find 

cCn\,-^ee^ = e^ (C.7) 

as expected from a naive estimate. 
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C.2 Misalignment from /5^ 

We perform the same analysis on the terms proportional to d, which implicitly encode the 
split between terms that carry factors of /f,^ versus /^^ from internal propagators. For Cy^ 
we have 

Crbld ~ {yi2 + Sy22) -7{yn +Sy2i). (C.8) 

Following the argument that the terms should cancel when the sign is swapped and using 
this to estimate the size of each bracketed term, one finds d Cji)\^ ~ e^, so that the net 
contribution of the d term is subdominant to (C.7). 

On the other hand, the /^^ misalignment in the C^^, term cannot be neglected, 

C'n\d ~ V {5'yi2 + y22) - (<J'yii +2/21) • (C.9) 

Here both terms are 0{e) so that the total contribution is 



dc: 



4 



3 

r^ t 

d 



~ e^ (CIO) 



nb 
which dominates over the term proportional to c in (C.5). 

C.3 Size of misalignment coefficients 

Thus the final order of magnitude estimate for the Cyf, and C^^ coefficients are 

Cn-cCnl^e^ (C.ll) 



75 ~ u '-^7(, 



Ov/i ^^ (i O' 



e% (C.12) 



so that unlike the anarchic contribution, the right-to-left (unprimed) and left-to-right 
(primed) Wilson coefficients are of the same order of magnitude. 

D Comments on 5D dipole theory uncertainties 

Finite 5D loop effects carry subtleties associated with cutoffs and UV sensitivity^. While 
the one loop contribution discussed in this paper is manifestly finite, higher loops are 
potentially divergent and require explicit calculations. Here we focus on the sensitivity of 
the finite loop-level result to UV physics at, for example, the strong coupling scale where the 
5D theory is expected to break down. In [37] it was pointed out that the naive dimensional 
analysis (NDA) for a brane and a bulk Higgs differ due to the dimension of the Yukawa 
coupling — the NDA two-loop contribution for the former gives an 0(1) correction relative 
to the one loop result, whereas this is not expected for the latter. In this appendix we 
comment on subtleties coming from 5D Lorentz invariance that may plausibly avoid this 
'worst case' NDA estimate. Indeed, the NDA for the one-loop contribution to these dipole 
operators is logarithmically divergent; one may understand the correct one-loop finiteness 
as coming from 5D Lorentz symmetry. 



^We thank K. Agashe, J. Hubisz, and G. Perez for discussions on these subtleties. 
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These comments are meant to demonstrate non-trivial points in these calculations that 
require particular care when drawing conclusions about UV sensitivity in these processes; 
a more careful investigation with explicit calculations of these effects is beyond the scope 
of this work. 

Note that the general features of the phenomenological picture presented in Section 6 
are unchanged even if there are 0(1) corrections to the Wilson coefficients. 

D.l KK decomposition 

5D Lorentz invariance imposes that in the KK reduced theory, the 4D loop momentum 
cutoff should be matched to the number of KK modes in the effective theory. This was 
mentioned in [39] to motivate a manifestly 5D calculation by pointing out that naively 
taking the finite 4D loop cutoff to infinity drops terms of the form (nMxK/A)^, where 
tiMkk is approximately the mass of the n KK mode. Indeed, from the 4D perspective 
this may appear to suggest a non-decoupling effect where the dominant contribution comes 
from heavy KK states so that the calculation seems to be sensitive to UV physics. 




Figure 13. A sketch of the 5D momentum space where the circle of radius A represents the 
boundary of a 5D Lorentz invariant loop momentum integration region. Marks on the kz axis show 
the masses of KK states. Dashed lines demonstrate that the 4D loop cutoff which respects 5D 
Lorentz invariance depends on the particular KK mode. 

However, as demonstrated in Fig. 13, imposing 5D Lorentz invariance requires that 
each KK mode carries a different 4D momentum cutoff. In particular, the n KK mode 
carries a smaller 4D cutoff A„ than that of the first KK mode, Ai since the momentum 
integral must fall within the circle of radius A, the 5D momentum space cutoff. Thus in 
4D the high KK modes are not sensitive to the same cutoff as lower KK modes. This gives 
a sense in which 4D decoupling can manifest itself while preserving 5D Lorentz invariance. 
In this sense it is difficult to use this matching to diagnose UV sensitivity. 

As a qualitative and demonstrative estimate, one can use the expression in Section 6.6 
of [39] for a neutral Higgs diagram and impose a KK number dependent cutoff for each 
state in the loop so that 5D Lorentz invariance is imposed as in Fig. 13. One finds that, 
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for example, in a sum of 200 KK modes, the highest 20 modes only contribute ~ 20% to 
the total result. 

D.2 5D cutoff 




A = 5 TeV 



y=R'k, 



Figure 14. Plot of the charged Higgs integrand as a function of the dimensionless loop momentum 
in the position/momentum space picture. The dashed line is a heuristic 5D cutoff A representing 
the strong coupling scale. The shaded region represents the error from taking the loop momentum 
to infinity rather than A; the contribution of this shaded region is approximately 15% of the total 
integral. 

Another way to diagnose UV sensitivity is to consider the effect of a cutoff in the 5D 
picture, for example, by setting a cutoff at A = 5 TeV representing the strong coupling scale 
at which the 5D theory breaks down. Fig. 14 shows the dimensionless integral associated 
with the charged Higgs loop, where y = R'kE is the dimensionless variable representing 
the loop momentum. Observe that the dominant contribution to the effect does not come 
from arbitrarily large y but rather in the peak at low values of y. Cutting off the integral 
at A = 5 TeV (dashed line) gives an error of approximately 15%, which is comparable to 
the subleading diagrams that were not included in this analysis. 
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